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Executive  Summary 
Introduction 

A  data  file  from  the  Experiment  on  Rapidly  Intensifying  Cyclones  over  the 
Atlantic  (ERICA)  was  obtained  and  graphically  rendered  on  a  Sun  workstation 
using  the  framework  of  a  particle  rendering  system.  The  mathematics  behind  the 
rendering  is  discussed  by  deriving  the  required  parametric  equations.  Techniques 
for  enhancing  the  graphical  display  of  the  cloud  liquid  water  content,  such  as 
dithering  and  jittering,  are  described  as  well.  The  pseudocode  for  rendering  is  also 
given,  illustrating  the  use  of  the  parametric  equations  and  the  enhancement 
techniques. 

Purpose 

The  goal  of  this  effort,  was  to  write  C  language  software  which  would  depict,  as 
realistically  as  possible,  a  cloud  image  from  cloud  liquid  water  content  data  on  a 
UNIX  workstation  such  as  a  Sun  or  a  Hewlett-Packard. 

Overview 

A  literature  search  was  first  carried  out  to  determine  what  work  had  been  done  in 
the  area  of  particle  systems.  Particle  systems  have  been  used  to  model  entities 
such  as  fire,  clouds,  and  water.  A  dramatic  example  of  their  use  is  in  the  movie. 
Star  Trek  II:  The  Wrath  of  Khan,  in  which  fire  and  explosions  are  very  realistically 
simulated.  [4]  Reeves  uses  stochastic  processes  to  generate  particles  in  a 
particular  volume  and  assigns  initial  velocities  to  each  particle.  These  concepts 
were  considered  until  the  author  found  that  clouds  could  be  rendered  using  simpler 
means  through  the  use  of  perspective  projection  equations  only.  [3] 

In  the  particle  rendering  system,  each  point  is  considered  as  a  light  source,  and 
where  the  point  projects  onto  the  screen  is  a  function  of  the  data  point  coordinates 
as  well  as  the  chosen  coordinates  of  the  eye.  Both  random  jitter  and  dithering  are 
applied  to  the  projected  points.  [3]  The  application  of  random  jitter  eliminates  the 
lattice  effect  one  sees  when  the  gridded  data  points  are  rendered  without  any 
massaging.  Dithering  provides  the  viewer  with  a  sense  of  magnitude  at  each 
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point.  This  is  accomplished  by  applying  a  dithering  matrix,  in  this  case  a  3  x  3 
matrix,  at  each  point;  whereby,  the  number  of  pixels  turned  on  within  that  matrix 
will  be  proportional  to  the  magnitude  of  the  data  value  at  that  point. 

Conclusions 

A  three-dimensional  rendering  of  cloud  liquid  water  content  can  be  done  fairly 
simply  on  a  UNIX  platform,  such  as  Sun,  with  the  appropriate  parametric 
equations  and  rendering  techniques  and  using  the  display  capability  of  the  XView 
library  of  routines  provided  by  Sun  Computer. 

Recommendations 

For  the  requirement  of  a  static  display  of  cloud  liquid  water  content  using  Sun’s 
XView  library  of  routines,  the  techniques  described  herein  should  prove  to  be 
quite  useful. 
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1.  Introduction 


1.1  Background 

A  technique  is  given  for  displaying  cloud  liquid  water  content  data  in  three 
dimensions.  The  data  file  from  the  Experiment  on  Rapidly  Intensifying  Cyclones 
over  the  Atlantic  (ERICA)  was  output  from  the  Limited  Area  and  Mesoscale 
Prediction  System  (LAMPS),  developed  at  Drexel  University,  which  produces  a 
three-dimensional  grid  of  the  original  data  set.  In  a  particle  rendering  system, 
each  point  within  a  three-dimensional  data  set  is  projected  onto  the  x-y  plane 
(computer  screen)  using  the  perspective  projection  equations.  The  rendering 
software  is  written  in  C  and  uses  the  XView  library  of  graphic  routines  from  Sun. 

1.2  Purpose 

The  purpose  of  this  report  is  to  describe  a  relatively  simple  technique  for  three- 
dimensional  display  of  cloud  liquid  water  content.  Using  Sun’s  XView  library  of 
routines,  this  technique  is  easily  implemented. 

As  this  was  an  independent  study  project,  the  author’s  advisor  clearly  stipulated 
that  commercial  rendering  software  packages  not  be  used,  but  rather  the  code 
should  be  “built  from  the  ground  up.”  The  effort  began  in  May  1994  when  the 
author  registered  for  an  independent  study  project  which  concludes  with  this 
documentation. 

1.3  Overview 

First,  the  ERICA  data  set  is  described  including  how  it  was  created  and  its 
dimensions.  Next,  the  techniques  used  in  a  particle  rendering  system  are  analyzed 
including  the  perspective  projection  equations,  as  well  as  dithering  and  jittering. 
The  Results  section  (2.3)  shows  the  graphical  display  of  cloud  liquid  water  content 
with  the  particular  parameters  used  to  achieve  the  image.  Lastly,  the  graphical 
display  software  is  outlined  in  terms  of  pseudocode. 
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2.  Cloud  Particle  Rendering  System 

2.1  Cloud  Liquid  Water  Content  Data 

2.1.2  ERICA  Data  Set 

The  data  set  used  for  this  project  is  an  83  x  53  x  34  (x,  y,  z)  grid  of  cloud  liquid 
water  content.  The  original  data  set  of  cloud  liquid  water  content  comes  from  the 
ERICA  experiment,  and  the  grid  is  output  from  the  atmospheric  dynamics 
software  module  called  the  LAMPS,  developed  at  Drexel  University.  The  areal 
extent  of  the  grid  is  from  99.658°  to  34.938°  West  longitude  and  from  22°  to 
54.5°  North  latitude.  The  layers  in  the  vertical  are  at  the  surface,  .025  km,  .375 
km,  .75  km,  then  every  0.5  km  up  to  15.25  km,  with  the  top  layer  at  16  km. 

2.2  Techniques 

2.2. 1  The  Perspective  Projection  Equations 

The  perspective  projection  equations  allow  one  to  project  points  from  x,  y,  z  space 
onto  the  computer  screen.  [1] 
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Figure  1.  Graphical  depiction  of  world  and  screen  coordinate  system,  eye 
position  (xe,ye,ze)  and  a  projected  point  (xp,yp). 
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The  equations  are  derived  as  follows: 


x-screen  =  x-axis  of  computer  screen, 
y-screen  =  y-axis  of  computer  screen. 

x-world  =  x-axis  for  liquid  water  content  data  set. 
y-world  =  y-axis  for  liquid  water  content  data  set. 
z-world  =  z-axis  for  liquid  water  content  data  set. 

d  =  distance  between  the  planes  formed  by  x-world,  y-world  and  the  screen. 

xi  =  x-coordinate  of  liquid  water  content  data  set. 
yj  =  y-coordinate  of  liquid  water  content  data  set. 
zk  =  z-coordinate  of  liquid  water  content  data  set. 

xe  =  x-coordinate  for  the  viewer's  eye. 
ye  =  y-coordinate  for  the  eye. 
ze  =  z-coordinate  for  the  eye. 

xp  =  x-coordinate  of  the  projected  point  on  the  computer  screen, 
yp  =  y-coordinate  of  the  projected  point, 
zp  =  z-coordinate  of  the  projected  point. 


X  =x  +p  *  (x.  -  x^) 

p  e  ^  ^  i  e  ' 

(1) 

1 

* 

+ 

II 

(2) 

(3) 

d=z^+p*{zj^-z^) 

(4) 
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The  equations  above  represent  the  parametric  equations  for  projecting  the  x,y,z 
point  in  the  liquid  water  content  data  set  onto  the  computer  screen. 

On  a  line  formed  by  the  points  (xe,  ye,  ze)  and  (xi,  yj,  zk)  p  represents  the 
fractional  distance  that  (xp,  yp)  is  from  (xe,  ye,  ze). 


As  you  can  see,  we  have  let  d  — 
At  d  =  0, 


z 


p 


p=- 


(5) 


so. 


-  z 


X  =x  + 
p  e 


Z.  -  z 

k  e 


*  (x.  -JC  ) 

V  I  e  ■' 


=  JC  + 
e 


X  *  Z  -  JC.  *z 

e  e  I  e 


z.-z 

k  e 


(6) 

(7) 


X 

e 


*Z,-JC  *Z  -JC.*Z  +x  *z^ 
k  e  e  I  e  e  e 


(8) 


X  *Z.-X.*Z 
e  k  I  e 


Z,  “Z 
k  e 


(9) 
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Similarly  in  the  y-direction, 


ye*^e-yj*^e 


(10) 


(11) 


(12) 


z.-z 

k  e 


(13) 


2.2.2  Jittering 

Random  jitter  is  applied  to  each  projected  point  by  assigning  one  of  the  eight 
immediate  neighboring  pixels  (or  xp,  yp)  to  that  point.  [3]  Its  purpose  is  to  avoid 
the  "lattice"  effect,  where  the  points  display  in  a  strict  grid-like  pattern,  unlike 
what  one  sees  in  nature.  This  is  accomplished  in  the  C  code  by  creating  two 
temporary  arrays  for  both  x  and  y,  each  of  size  3: 

yptemp[0]=yp; 
yptemp[l  ]=yp-spacing; 
yptemp[2]=yp+spacing; 

xptemp[0]=xp; 
xptemp[  1  ]=xp-spacing; 
xptemp[2  ]=xp+spacing; 


14 


Spacing  is  a  random  displacement  from  the  centeipoint  and  can  be  set  by  the  user. 
The  actual  point  to  be  rendered  will  be  selected  by  creating  a  value  randomly  of 
either  0, 1,  or  2  by  using  drand48(),  a  built-in  C  function.  The  resulting  value  is 
used  as  an  index  into  the  arrays  described  above.  The  point  rendered  can  be  the 
original  point  that  was  projected  or  any  one  of  eight  random  displacements  around 
that  point. 


2.2.3  Dithering 


Dithering  provides  a  means  of  signifying  magnitude.  A  3  x  3  dithering  matrix  is 
used  which  provides  for  ten  levels  of  intensity.  [2]  The  values  of  the  liquid  water 
content  range  from  0  to  0.83  g/cml  This  range  was  split  into  10  intervals.  Thus, 
the  value  of  the  concentration  at  a  pixel  determines  how  many  pixels  are  turned 
on  within  that  3x3  neighborhood.  The  patterns  adopted  in  this  case  are  as 
follows: 


liquid  water  content  range  pixel  pattern 

1.  0.0-  .08  [  ] 

[  ] 

[  ] 

2.  .09 -.16  [  ] 

[  *  ] 

[  ] 

3.  .17 -.25  [  ] 

[*  *  ] 

[  ] 

4.  .26 -.33  [  ] 

*  ] 

[  *  ] 
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5.  .34 -.42 


[  ] 

j-*  * 

[  *  ] 


6.  .43 -.50 


7.  .51 -.59 


8.  .60 -.67 


9.  .68 -.76 


10.  .77 -.83 


[  *] 
[*  *  *] 
[  *  ] 

[  *] 

r*  *  *1 

[**  ] 

[**  ] 

[*  *] 
[*  *  *] 
p  * 

[*  *  *] 

|^=C  >«:  4=J 
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2.3  Results 


Figure  2.  Three-dimensional  display  of  ERICA  liquid  water  content. 


Displayed  above  is  the  output  of  the  graphical  rendering  of  the  ERICA  liquid 
water  content  data  set  as  preprocessed  by  the  LAMPS  model.  The  lowest  grid 
layer  is  easily  distinguishable  in  the  lower  right  hand  portion  of  the  figure  slanting 
from  southwest  to  northeast. 

Each  point  in  the  cloud  (xi,yj,zk)  is  treated  as  a  light  source.  Its  corresponding 
perspective  projection  point,  (xp,yp),  is  the  point  that  is  illuminated.  In  this 
particular  case,  xi  ranges  from  1  to  83,  yj  from  1  to  53,  zk  from  1  to  34.  It  is 
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assumed  that  the  left  eye  and  the  right  eye  are  at  the  same  altitude  and  latitude  and 
are  separated  by  an  interocular  distance,  e. 


That  is: 


^  T7T  ^  IP  t>  ^ 

EL  ER  e 


(14) 


y  EL  y  ER  y  t 


(15) 


^ER  ^EL  ^ 


(16) 


where  EL  refers  to  the  left  eye  and  ER  to  the  right  eye. 


The  disparity,  d,  is  defined  to  be: 


d  Xpp 


(17) 


and  represents  the  difference  in  xp  between  the  left  and  right  eyes. 


X 


e 


e 


have  been  set  to  15,  -2100,  12,  and  300,  respectively,  for  this  particular  rendering. 
These  first  three  values,  respectively,  place  the  viewer  near  the  western  edge,  back 
fi-om  the  cloud  looking  northward,  and  looking  down. 

A  threshold  for  the  minimum  value  of  all  the  liquid  water  content  array  values  to 
be  displayed  was  set  to  0.1 1. 


XSCALE,  YSCALE,  and  ZSCALE  are  parameters  used  when  the  current  grid 
point  is  being  projected.  Their  purpose  is  to  scale  the  data  so  that  it  will 
reasonably  fill  the  particular  frame  size  chosen.  The  values  used  were  5,5,  and 
-0.3  respectively.  These  positive  x  and  y  values  have  the  effect  of  expanding  the 
cloud.  See  section  2.4  for  the  particular  equations. 

XSHIFT,  YSHIFT,  and  ZSHIFT  provide  translations  in  the  x,  y,  and  z  directions. 
The  respective  values  used  were  400,  1400,  and  10.  Section  2.4  illustrates  how 
they  are  used  in  the  equations. 

The  data  file  from  the  ERICA  experiment  contains  149,564  data  points  total  (two 
points  serve  as  bookends)  with  12,368  non-zero  liquid  water  content  values.  The 
minimum  data  set  value  is  0.0,  and  the  maximum  value  is  0.83  g/cm^  of  liquid 
water. 

2.4  Outline  of  the  C  Source  Code  for  Rendering  Cloud  Liquid 
Water  Content 

2.4.1  function  mainQ 

READ  the  83  x  53  x  34  grid  of  liquid  water  content  (Iwc)  values  into  a  3-d  array. 
CREATE  a  labeled  window  frame  using  xv_create(). 

CREATE  colormap  segment  with  xv_create(). 

SET  frame  height,  width  with  a  call  to  xv_set(). 

CREATE  a  conversion  table  for  pixel  values  from  Xview  to  XI 1 . 

CREATE  a  canvas  within  the  frame  and  render  the  Iwc  values  by  calling 
xv  createO  which  in  turn  calls  canvas_repaint_proc(),  the  procedure  which  does 
the  rendering. 
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2.4.2 


procedure  canvas_repaint _procQ 


LOOP  over  the  34  levels  (k) 

SCALE  the  current  z  value  for  the  frame  size  as 

=  ZSCALE  *k+ZSHIFT  (18) 


DETERMINE  the  disparity  as 


7.  ^  P. 


LOOP  over  the  53  columns  (j) 

SCALE  the  current  y  value  for  the  frame  size  as 

y.  =  YSCALE  *7  +  YSHIFT  (20) 


DETERMINE  where  the  y  point  will  project  for  rendering  as 


Yr 


(Ye  *  -  (Yj  * 


(21) 


LOOP  over  the  83  rows  (i) 

SCALE  the  current  x  value  as 


X .  =  XSCALE  *  i  +  XSHIFT 

I 


(22) 


20 


DETERMINE  where  x  point  will  project  for  left  eye  as 


(a:,  *Zt)  -  (x, 


(23) 


DETERMINE  where  x  point  will  project  for  right  eye  as 
^PR  ^  )Xpj^-  d  (24) 

JITTER  the  current  data  point. 

DRAW  the  appropriate  dither  pattern  given  the  value  of  the  current 
data  point. 

END  LOOP 

END  LOOP 

END  LOOP 
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3.  Conclusions 


A  fairly  realistic  static  depiction  of  a  cloud  was  produced  from  the  cloud  liquid 
water  content  data,  A  time-lapse  sequence  depiction  of  a  cloud  would  capture  its 
dynamic  nature,  and  this  static  depiction  is  a  good  first  step  toward  that  goal. 
Further  experimentation  with  different  dithering  matrices  and  different 
perspectives  would  be  beneficial  to  see  if  any  improvements  could  be  made  to  the 
system. 
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4.  Recommendations 


A  conversion  from  Xview  to  Motif  should  be  carried  out,  since  Motif  is  the 
government  standard  used.  It  is  also  unclear  if  Sun  will  continue  to  support 
Xview. 

The  author  has  had  recent  discussions  with  Professor  McCoy  on  ways  to  make  the 
cloud  rendering  more  realistic.  One  method  discussed,  would  render  cubes  rather 
than  just  pixels.  The  color  for  a  particular  cube  face  would  be  a  funetion  of  the 
amount  of  water  droplets  between  the  reflecting  cube  and  the  eye,  and  the 
magnitude  of  the  liquid  water  content  at  that  point.  A  light  source  could  also  be 
added,  which  would  involve  taking  shadows  into  account. 
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